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ABSTRACT 

The strong time-dependence of the dynamics of galactic bars yields a complex and rapidly 
evolving distribution of dense gas and star forming regions. Although bars mainly host re¬ 
gions void of any star formation activity, their extremities can gather the physical conditions 
for the formation of molecular complexes and mini-starbursts. Using a sub-parsec resolu¬ 
tion hydrodynamical simulation of a Milky Way-like galaxy, we probe these conditions to 
explore how and where bar (hydro-)dynamics favours the formation or destruction of molec¬ 
ular clouds and stars. The interplay between the kpc-scale dynamics (gas flows, shear) and 
the parsec-scale (turbulence) is key to this problem. We find a strong dichotomy between the 
leading and trailing sides of the bar, in term of cloud fragmentation and in the age distribution 
of the young stars. After orbiting along the bar edge, these young structures slow down at the 
extremities of the bar, where orbital crowding increases the probability of cloud-cloud colli¬ 
sion. We find that such events increase the Mach number of the cloud, leading to an enhanced 
star formation efficiency and finally the formation of massive stellar associations, in a fashion 
similar to galaxy-galaxy interactions. We highlight the role of bar dynamics in decoupling 
young stars from the clouds in which they form, and discuss the implications on the injection 
of feedback into the interstellar medium, in particular in the context of galaxy formation. 
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1 INTRODUCTION 

In isolated disc galaxies in the local Universe in general and the 
Milky Way in particular, most of the star formation activity oc¬ 
curs in giant molecular clouds distributed within kpc-scale struc¬ 
tures like spiral arms, where the interstellar medium (I SM) is 
about one hundred time d enser than in inter-arm regions l lVisseil 
Il980l : lllughes et al.ll201.'lh . Such excess of star formation is also 
observed in the central regions of galacti c bars, again bec ause 
of high surface densities o f molecular gas dsheth et al.ll2002L but 
see lUongmore et al.l 1201.31 on a deficit of star formation in the 
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central molecular zone). Furthermore, the extremities of bars can 
host large amounts of molecular gas, distributed i n massive com¬ 
plexes, also called giant molecular associations l|Solomon et al.l 
Il987l : l^artin & FriedlilllQQTl : iNguven Luong et al.ll201 1 1. Among 
them, the complex W43 at the tip of the Milky Way bar gath¬ 
ers about 9 X 10® M© of ga s, embedded in a neutral hydro¬ 
gen envelope of 3 x 10® Mq jMotte et al.ll2014l) . This complex 
hosts a mi ni-starburst activity, with the ongoing f ormation of mas¬ 
sive stars jMotteetalj 20031 : Louvet et aiTl2014h . Using position- 
velocity data, Motte et ^ ( l2014h proposed that such activity results 
from a collision between clouds accreted from a nearby spiral arm. 
However, because of the very short dynamical timescales in such 
environment (~ 10 Myr) and the interplay between the numerous 
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physical processes at stake, understanding the origin of molecular 
complexes, their lifetimes, and the physical conditions leading to 
enhanced star formation in bars is rather difficult observationally. 

Numerical simulations provide a complementary approach 
to the problem, and have been long used to address the ques¬ 
tions related to the dynamics a nd star formation in bars (see 
a rev iew in lAthana ssoulal l2 013l) . Using sticky particle simula¬ 
tions, ICombes & GerinI i 19851) noted asymmetries in the distribu¬ 
tion of dust lanes between the leading and trailing sides of bars. 
Suc h results have also been established in hydrodynamica l mod¬ 
els ) |Prendergasjl983l : lAthanassoulall992l : |Piner et kjl995h which 
also showed that gas concentration on the leading side corresponds 
to shocks. These simulations have sufficiently high resolutions 
(~ 100 pc) to describe the inner structures of the bar. However, 
such pioneer studies faced technical limitations and could not re¬ 
solve how these features fragment (or not) into molecular clouds 
(~ 10 pc). Although modem simulations still face the long-lasting 
issue of coupling between the kpc-scale dynamics of the galaxy and 
the parsec and sub-parsec physics driving star formation and the as¬ 
sociated processes like feedback, technical improvements in high 
performance computation push further the limitations. Present-day 
galaxy simulations now encompass both scales at once (s ee e.g. 


galaxy simulations now encompass bom scales at once (s ei 

IPobbs et al.ll2012l: iBonnell et al.ll2013l: iFuiimoto et al.ll2014 f 


_ Using the model of the Milky Way presented in Renaud et al.l 

ll2013h . we address here some of the questions related to the evo¬ 
lution of the structure of the ISM and the formation of stars in a 
galactic bar, focussing on the regions at their tips. The main advan¬ 
tage of this simulation is to describe the evolution of the interstellar 
medium and the formation of stars by resolving the giant molecu¬ 
lar clouds at parsec and sub-parsec scales, within a “live” galactic 
context, i.e. with no forcing of the potential nor of the spiral or bar 
pattern. Therefore, this simulation encompasses both dynamical as¬ 
pects from the kpc scales like tides, shear, disc rotation, and local 
parsec and sub-parsec effects like cloud self-gravity, turbulence and 
stellar feedback. This paper aims at providing comprehensive sce¬ 
narios for the formation of gaseous structures, their fragmentation 
into molecular clouds, the interaction of these clouds, and the po¬ 
tential formation of molecular complexes like W43. 


2 NUMERICAL SIMULATION 

Details of the galax y simulation used in this work are presented in 
iRenaud et al.l 1 1201^ and summarised here. We haye run a hydrody¬ 
namical simulation of a Milky Way like gala xy in isolation, using 
the adaptive mesh refinement code RAMSES )lTevssieill2002h . The 
initial conditions are axisymmetric distributions of stars (bulge and 
discs) and dark matter halo rendered with particles, and gas mod¬ 
elled by an exponential disc on the mesh of the simulation. The 
stellar component for ms spirals and a bar within a few 10® yr (see 
lEmsellem et al.l2015h . The gas then naturally follows this potential 
well, and further fragments to form its own structures. The forma¬ 
tion and evolution of stellar and gaseous structures naturally arise 
from the dynamical evolution from axisymmetric initial conditions 
and are not regulated by a pattern speed nor a fixed potential. 

In a cubic box 100 kpc wide, the maximum resolution reached 
is 0.05 pc for a cloud crossing time (i.e. a few 10® yr), and sub¬ 
parsec resolution for a few cloud lifetimes (i.e. several 10 Myr). 
The equation of state, which is a fit to the actual thermal heat¬ 



Figure 1. Map of the stellar component of our Milky Way model (seen face- 
on) at the beginning and the end of our period of interest. (The bar rotates 
anti-clockwise in all our Figures.) 


ing/cooling balance, mimics a hot diffuse halo (10®“^ K), an 
isothermal warm phase (10'* K up to « 10“* cm“®) and a poly¬ 
tropic branch at high density, down to 14 K. In Milky Way-like 
galaxies, only a moderate amount of gas lies off the adopted equa¬ 
tion, and the heating and cooling processes make it rapidly fall 
back onto the relation, such that this equation of state provides a 
go od approximation to more elaborate implementations, as shown 
in iKraliic et alJ bOldh . A pressure floor (for gas denser than « 
2 X 10® cm“®) ensures that the local Jeans length is always re¬ 
solved by a few cells. 

Star formation proceeds at a constant efficiency per free-fall 
time in the gas denser than 2000 cm“® (which corresponds, on 
average over the galaxy, to the transition where self-gravity takes 
over the supersonic turbulent dominated regime, visible as a power- 
law tail in th e otherwise log-norm al density probability distribution 
function, see 'Renaud et alj2013L their Figure 10). The prescription 
for stellar feedback consists in various phenomenological recipes 
for photo-ionisation in Hll regions (by setting the temperature of 
Stromgrem spheres to 2 x 10* K), radiative pressure (by deposi¬ 
tion of radial velocity kick s in these spheres) and kinetic supernova 
blasts jRenaud et al.l201^ . No magnetic field has been included. 

The snapshot at f = 800 Myr (with t = 0 corresponding to 
the initial conditions of the simulation) provides a good match to 
the observational data of the Milky Way, in term of morphology and 
kinematics. Although some fine features are absent of the model 
(e.g. the nuclear ring at 200 pc), the main str uctures of the Milky 
Way (bar, spiral arms) are well reproduced (see lRenaud et al.l2013l : 
lEmsellem et al.llToi5l : IChemin et alj|2015l) . In that sense, although 
it is still utopian to compare the simulation with observations at 
the level of individual GMCs, this model can confidently be used 
to probe the average physical conditions of the real Galaxy and 
address the questions of formation and evolution of GMCs from a 
statistical perspective. 

In this contribution, we study the formation of gas clouds and 
stars in the outer parts of the bar. The (hydro-)dynamics of the in¬ 
nermost kpc, and the coupling b etween the large- s cales a nd the cen¬ 
tral black hole are presented in lEmsellem et alJ l l2015l) . We focus 
our analysis between t = 755 Myr and t = 800 Myr to monitor 
the physical conditions leading to the formation of clouds and then 
to stars (Fig. [T). Running the simulation at sub-parsec resolution 
over this period would be too costly. Instead, the resolution is pro¬ 
gressively increased, from 3 pc to 0.05 pc, one refinement level at 
a time. Before activating a finer level, we let the simulation run for 
several free-fall times of the densest structures, to ensure the signa- 
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Figure 2. Map of the gas density (at t = 755 Myr), in the bar and the inner 
region of the disc. 


tures of previous refinements have been erased. Star formation and 
feedback are however active throughout the analysed period. Note 
that the effective resolution controls the fragmentation level of the 
ISM, which enables us to allow for the formation of spirals and bar 
while preventing that of molecular clouds during the earliest epochs 
of the simulation 3-4 rotation periods) to avoid a too floccu- 
lent disc. The coarsest resolution is however sufficient to describe 
the properties of the ISM at the scale of molecular clouds (several 
10 pc), but without being able to describe their inner structures. 

The bar formed at f « 600 Myr and evolved from a weak 
bar (maximum amplitude of the tangential to radial force ratio 
of Q « 0.2 at t = 700 Myr) to a stronger bar (Q « 0.6 at 
t = 800 Myr, see lEmsellem et al.]|2015l . their Section 3.1). At 
t = 800 Myr, the bar pattern-speed is « 58km/s/kpc (with 
co-rotation found at 3.6 kpc). Its semi-major axis is 2.4 kpc, with 
an axis ratio of 1:3. Inner Lindblad resonances, as computed from 
the azimuthally averaged mass distribution, are found at 40 pc and 
450 pc from the galactic centre. More details on the evol ution 
of the disc and the ba r can be found in iRenaud et al.l ( 1201 3h and 
lEmsellem et al.l ( l2015h . 


Figure 3. Map of the turbulent Mach number evaluated at the scale 24 pc. 
Dashed lines indicate the major and minor axes of the bar, and thus tell 
apart the leading and trailing quadrants. The highest turbulent Mach num¬ 
bers are found on leading edges of the bar, where the ISM gets supersoni¬ 
cally shocked. 


leading trailing leading trailing 



Azimuth [°] 


Figure 4. Azimuthal profile of the turbulent Mach number (computed at the 
scale of 24 pc), along the edge of the bar. The edge is here defined as an el¬ 
lipse of semi-major and semi-minor axes 1.8 kpc and 0.7 kpc respectively, 
and with an inclination of 72 degrees with respect to the axis x = 0 (see 
Fig.H). To account for the non-perfectly ellipsoidal nature of the bar and 
include the densest gas features, data is extracted in a ring spanning 250 pc 
on both the inner and outer sides of this ellipse. Azimuths are indicated on 
Fig.H] On each quadrant (matching those displayed in Fig. [3), a horizontal 
dashed line indicates the average turbulent Mach number. The leading sides 
of the bar host regions with higher turbulent Mach numbers than the trailing 
quadrants. 


3 RESULTS 

Fig.ll shows the density of gas at f = 755 Myr, before most of 
the star forming clumps form in the bar. Using this snapshot al¬ 
lows us to measure the physical conditions of cloud formation. Gas 
has already gathered in dense structures, following the gravitational 
potential of the stellar component, and is at the turning point of 
fragmenting into clouds. Some regions along the outermost spiral 
arms (e .g. on lower part of Fig. |3 already host beads-on- a-string 
clouds dElmegreen & Elmegreenll983l : lRenaud et aljTtlTsI) . which 
are still in the phase of accreting their surrounding material. Gravi¬ 
tational torques from the bar and resonances regulate the gas flows 
toward the centre a long mini-spirals, in pa rticular inside the inner 
kpc, as discussed in lEmsellem et al.l ( l2015l) . 


3.1 Tlirbulence 

At the scale of 24 pc, only a few regions yield supersonic tur- 
bulence (Fig. [3]l : the leading edges of the bar (as observed by 
ISheth et al.ll200il and, to a lower extent, the two main spiral arms 
close to the points where they c onnect to the bar. ( The turbulent 
Mach number is computed as in lKraliic et al.ll20l3 . their Section 
4.1.) To further illustrate this point. Fig. |4] displays the azimuthal 
profile of the turbulent Mach number along the edge of the bar, 
confirming the presence of a leading versus trailing asymmetry. 

These high turbulent Mach numbers mostly originate from a 
velocity dispersion (~ 20 km s“^) about 2 to 10 times larger in 
these regions than elsewhere along the bar edge (where density, 
temperature, and thus sound speed are comparable). The locations 
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Figure 5. Ratio of the compressive turbulence energy (encompassing both 
compression and rarefaction) to the solenoidal turbulence energy, both eval¬ 
uated at the scale of 24 pc. The dotted ellipse and lines indicate the approx¬ 
imate location of the densest gas structures along the edge of the bar and the 
spiral arms. The horizontal mark in the color bar indicates energy equiparti- 
tion between the two turbulent modes (i.e. a ratio of 1/2, see text for details). 


of these regions suggest the presence of a supersonic shock be¬ 
tween the edge of the rotating bar 60 km/s/kpc, 50-500 K, 
~ cm“®) and the warm neutral medium (2500-3000 K, 

~ 10 cm”®) surrounding it. (Note that hotter and more diffuse 
inter-arm medium is found at larger radii, and only faces the spi¬ 
rals and not the bar itself.) Oppositely, the inner region of the bar 
stands in transonic and subsonic regimes. 

Number of studies revealed that, on top of the turbulent Mach 
number which describes the level of turbulence of the ISM, the 
nature of the turbulence can also be key in the formation and evolu¬ 
tion of dense gas structures (see e.g. lWada et alj2002l: Kritsuket_^ 
l2007l : [^derrath et alj20ld at cloud-scales, and Renaud et al.ll2014 
at galactic scales). This complementary aspect can be probed by 
measuring the relative importance of the compressive (curl-free) 
and solenoidal (divergence-fre e) modes of turbulen ce. To do so, we 
follow the method presented in iRenaud et al.l l l20l4 their Appendix 
A). The divergence and curl of the velocity field are computed with 
first-order finite differences using a group of 2® grid cells. In such 
scheme, the large scale component of the total velocity field (i.e. 
the non-turbulent term) cancels out, and only the turbulent compo¬ 
nents of the velocity field remain. By multiplying the divergence 
(respectively curl) by the scale-length over which it is evaluated 
(i.e. the cell size), one reconstructs the compressive (respectively 
solenoidal) velocity, which then can be expressed as a kinetic en¬ 
ergy. Note that the curl-free component represents both a compres¬ 
sion and a rarefaction effect, but we refer to it as compressive mode, 
for short. 

Fig. a maps the ratio of the energies of the two turbulent 
modes. Because the compressive mode carries one degree of free¬ 
dom while the solenoidal mode has two, energy equipartition be¬ 
tween the two modes is reached when their ratio equals one half. 
The discrete nature of the grid on which the motions are repre¬ 
sented implies that this ratio is over-estimated if ev aluated over too 
few resolution elements (see lFederrath et apTOllI) . By varying the 


resolution at which we compute the map of Fig.|5] we ensured that 
convergence is reached at 24 pc, i.e. over at least 8 cellfl 

Sharp density gradients, across the spiral arms and at the edge 
of the bar, are often associated with a change from compressive- 
to solenoidal-dominated turbulence. Most of the dense structures 
(bar and spirals) yield turbulence close to equipartition, or slightly 
solenoidal dominated, indicating the presence of local shear mo¬ 
tions. Only the low density gas (mostly in the inter-arm regions) 
yield a compressive-dominated turbulence (rarefaction) over kpc- 
scale volumes. Recall that at this stage of the simulation, star for¬ 
mation has not taken place yet and thus, this turbulence is not gen¬ 
erated by feedback, but only by the large-scale (hydro-)dynamics. 
Later, the turbulence pumped by stellar feedback at parsec scales 
will combine with that injected at kpc-scale by the galactic dynam¬ 
ics. 

Along the bar edge, the leading versus trailing dichotomy 
noted in Fig. and |4] for the Mach number is not recovered in the 
nature of turbulence. The regions along the edge share comparable 
turbulence as the spiral arms. No major differences are found be¬ 
tween the edge of the bar and the central kpc. The main differences 
in the organisation of the ISM is thus originating from large-scale 
aspects like shear (see next Section). 

Therefore, at the scale of molecular clouds, the overall turbu¬ 
lence intensity (Mach number) shows a rather strong dependence 
on the galactic environment, leading to differences in the struc¬ 
ture of the ISM and the further formation of dense clouds. The na¬ 
ture of turbulence (compressive versus solenoidal) in the dense gas 
structures is however fairly independent of the kpc-scale structures, 
contrarily to interac ting galaxies where it is mostly compressive 
llRenaud et al.l2014h . The reason for this might be that compressive 
tides, which are thought to drive large-scale compressive turbulence 
in interacting galaxies, do not exist over kpc-scale volumes in iso¬ 
lated galaxies. 

This picture remains valid at later stages in the simulation, 
even after the formation and collapse of clouds. 


3.2 Large-scale velocity field 

Orbi tal motions inside the bar induce a shearing effect on the ISM 
(see lAthanassoulalll992l : lEmsellem et aLllroT^ . Such shear pre¬ 
vents the formation of dense structures like clouds, which later 
translates in a lack of star formation inside the bar, despite rela¬ 
tively high densities. At the edge of the bar however, supersonic 
turbulence (as discussed above) balances a weaker shear (and tides) 
than in the innermost regions, and favours the formation and sur¬ 
vival of overdensities. 

Fig.[6]shows the maps of the radial and tangential components 
of the gas velocity in the reference frame of the bar. The highest 


^ In the lowest density regions (e.g. the inter-arm), the adaptive nature of 
the grid refinement implies that the cell size is larger than in denser zones 
like the bar and the spiral arms. Because of this, at fixed scale, the tur¬ 
bulent modes are computed over fewer resolution elements, and thus the 
importance of the compressive mode there (as a rarefaction effect) is over¬ 
estimated. In our case, only the extreme values in over-compressive regions 
(of very low density) are significantly changed (up to a factor of 3 in the 
most extreme cases), but the overall behaviour remains unaffected if quan¬ 
tities are computed at 100 pc scale instead of 24 pc. 
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Figure 6. Top: gas radial velocity. Blue and red indicate inward and out¬ 
ward motions respectively. Bottom: gas tangential velocity, minus the pat¬ 
tern speed of the bar (58 km/s/kpc). Blue and red indicate clockwise and 
anti-clockwise motions respectively. The dotted ellipse and lines indicate 
the approximate location of the densest gas structures along the edge of the 
bar and the spiral arms. 


radial velocities are found along the edge of the bar (pointing in¬ 
ward on the leading sides and outward on the trailing sides). Shear 
is clearly visible perpendicular to the bar major and minor axes 
(in particular in radial velocity) whereas kpc-size regions along the 
edge show almost uniform velocities. In (relative) absence of large- 
scale shear in the latter areas, small-scale (< 100 pc) motions dom¬ 
inate the velocity field, which translates into the turbulence patterns 
discussed in the previous Section. 

Fig.0displays the velocity component parallel to the bar mi¬ 
nor ax is along three slits, following the approach of lAthanassoulal 
1 I 1992 I her Figure 11). In the outermost regions considered here, 
over-densities correspond to spiral arms. Any of the associated ve¬ 
locity profiles within the arm (i.e. over ~ 100 pc) is monotonically 
decreasing, which indicates converging flow toward the arm centre. 
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Figure 7. Top: gas density map of the bar (as in Fig. [3- The black cir¬ 
cles show the inner Lindblad resonances (40 pc and 450 pc). The three 
dashed lines indicate the position of “slits” parallel to the minor axis, and 
shifted along the vertical axis by 500 pc, 0 pc and —500 pc. The bottom 
three panels show the velocity component parallel to the minor axis, along 
the three slits (with the convention of positive velocity representing motion 
from left to right). To ease the interpretation, the background greyscale im¬ 
age represents the logarithm of the gas density and the curve colour codes 
the local velocity gradient. With the sign convention adopted, red indicates 
a converging flow, blue a diverging flow, white and symmetric features de¬ 
notes bulk motions. A density map of the corresponding region is shown 
above each panel to help the identification of structures. 
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Figure 8. Profile of the velocity component perpendicular to the minor axis, 
along the three slits shown in Fig.|2] Background indicates the gas density 
in the slits, as in Fig.|2] 

At larger scale, the flows are set by the disc dynamics and thus vary 
with location. 

The innermost dense gas structures (only visible in the central 
slit) match the positions of the inner Lindblad resonances, at 40 pc 
and 450 pc. The gas is fuelled there along mini-spiral arms by the 
gra vitational torques of the bar and piled up o n these stable orbits 
(see lRenaud et alj2013l : lEmsellem et alj2015ll . Regions in between 
these orbits yield either low velocities or diverging flows, which 
participate in the lack of cloud and star formation in this area. 

At the edge of the bar (0.5 kpc < I®! 1 kpc), the densest 

structures match steep velocity gradients, i.e. shocks, especially on 
the leading side of the bar (a; « 300 ± 100 pc fo r SI and x ~ 
600 zb 100 pc for S3), which confirms the findings of lAthanassoulal 
( Il992h . On the trailing side, the velocity profile is mostly symmetric 
over a few 100 pc, which simply denotes the rotation of the bar. 

However, structures at intermediate radii and on the trailing 
side are not systematically associated with shocks. The difference 
originates from the velocity component perpendicular to the minor 
axis of the bar, of which profile is displayed on Fig. [8] Dense gas 
structures associated with shocks (i.e. with steep parallel velocity 
profiles on Fig.|7} are found in zones of low shear (i.e. flat perpen¬ 
dicular velocity profiles Fig.[^, particularly close to the edge of the 
bar, and in a few ~ 100 pc-scale regions where the perpendicular 
velocity profile locally flattens. 

Such relatively weak shear, combined with high densities 
and strong turbulence (recall Section IrTT l. shapes elongated dense 
gaseous structures and allows for their survival until they fragment 
into clouds. Most of the material found there does not spiral inside 
the bar but remains o n the bar edge and circulates around the bar 
(on x i-type orbits, see lContopoulos & Grosbolll989l : ISkokos et ^ 
l2002h at high velocity (up to 230 km s ^ in the bar reference 
frame). It takes about 20 Myr to travel from one of its apocen- 
tre to the next (on the other side of the bar, « 3 kpc on aver¬ 
age). The eccentric nature of this motion makes the gas slow down 
(~ 40 km s“^) when it reaches the apocentre, i.e. the extremities 
of the bar (along the major axis). Here, the spread in orbital en¬ 
ergies translates into a radial spread of apocentres, such that the 
material coming from the thin filaments parallel to the bar spans 
a large area (about 1 kpc^, between 1 kpc and 2 kpc away 



Figure 9. Age of the young stars in the vicinity of the tip of the bar at 
t = 800 Myr (i.e. 45 Myr after the instant showed in Fig. |^. Only 

one side of the bar is shown, for the sake of clarity. An arrow indicates the 
rotation of the bar. The object circled is discussed in Section lT^ The bottom 
panel shows a zoom-in of the rectangle area in the top panel. We consider 
three age bins: < 20 Myr (blue), between 20 and 40 Myr (green), and 
> 40 Myr (red). 

from the galactic centre). Despite this large volume in space, the 
slowing down of gas clouds makes them accumulate there (before 
being accelerated again along the oposite bar edge). We identified 
about 50 dense clouds (> 2000 cm“®) at each extremity of the 
bar, with a velocity dispersion between them of ~ 40 kms“^. 
This accumulation in space and time, so-called orbital crowding 
dKennev & LordflOOlh . increas es the probability of cloud-cloud 
collisi ons, as already noted by [Rodriguez-Femandez & CombesI 
(l2008h . Such events lead to star formation, as discussed in Sec¬ 
tion [T4] 


3.3 Young stars 

The natural consequence of the formation of dense clouds is that 
self-gravity will eventually take over regulating mechanisms and 
will lead to cloud collapse and star formation. 

Fig.m shows the spatial distribution and the age of the new 
stellar particles (i.e. the stars formed out ot gas during the simula- 
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Figure 10. Velocity field of the young stars, in the reference frame of the 
bar (58 km/s/kpc), in the same region and at the same epoch as in Fig. [9] 
The background displays the gas density map. Like the clouds, the stars 
slow down when reaching the exttemity of the bar. 



tion), about 45 Myr after the instant previously discussed, once the 
clouds have formed, collapsed and began to form their stellar ma¬ 
terial. The lack of star formation in the inner bar (except at the very 
centr e) is caused by the str ong shear destroying clouds in this re¬ 
gion temsellem et alj2015L see also Fig. D, and the orbital pattern 
preventing stars formed elsewhere to enter this area. 

One can tell apart three stellar populations: (i) the very young 
stars (< 20 Myr), (ii) the intermediate class (between 20 and 
40 Myr old), and (iii) the older population (> 40 Myr). (These 
classes are defined arbitrarily, only for the sake of clarity.) 

As expected, the very young stars are found associated with 
dense gas structures: along the spiral arm connected to the tip of 
the bar (i.e. in the elongated blue feature on the left-hand side of the 
figure), and in central regions of some dense clumps at the tip of the 
bar itself (bottom panel), with the notable exception of the massive 
clump on the edge of the bar [x « —0.1 kpc, y « —0.8 kpc) 
which will be discussed in Section [T4l 

The intermediate age stars (20 Myr < age < 40 Myr) are 
preferentially found at the tip of the bar, but on the trailing side 
(upper half of the bottom panel of Fig. H. Their ages indicate that 
they formed at the opposite tip of the bar and have travelled along 
the (upper) edge of the bar within the last ~ 20-40 Myr, at the 
high velocity noted above. When reaching this extremity of the bar 
(i.e. their orbital apocentre), these stars slow down, as shown in 
Fig.IIi They stay in the area of the tip of the bar for ~ 20 Myr 
before being slingshot along the leading side of the bar (bottom on 
Fig- m and [To) toward the other tip. This explains that stars on the 
leading side of the tip are a couple of 10 Myr older than those on 
the trailing side: they aged in the area of the tip. 

By tracking back the origins of the stars found at the extremity 
of the bar (by selecting all stars in the bottom panel of Fig. |9), we 
find that about 50% of them form in spiral arms and are accreted 
later by the bar, with the furthest formation site being 4.5 kpc away 
from the galactic centre. Among this accreted population, about 
half (55%) originates from the main spiral arms connected to the 
bar, while the other half forms along a secondary structure, equiv¬ 
alent to the 3-kpc arm (or Norma arm) in the real Galaxy. While 
these stars form in clusters, the dynamics of the accretion (shear 
and tides) makes most of these objects dissolve before they reach 
the extremity of the bar. They participate in populating the “field”. 



Figure 11. Gas density map in the region of interacting clouds along the 
edge of the bar (at t = 800 Myr). Signatures of cloud-cloud interaction are 
visible, as tidal tails and bridges made of dense gas. 


as opposed to the clustered distribution of the young stars formed 
in situ, as visible on Fig.j^ 

To summarise, the highly dynamical environment of the galac¬ 
tic bar directly influences the spatial distribution of young stars. It 
is however likely that the stellar populations will mix over several 
orbital periods (~ 100 Myr) and that the long-term evolution of the 
bar and the associated orbits would amplify this mixing process. 


3.4 Cloud-cloud collisions and W43-analogue 

The cluster of young stars at the position x « —0.1 kpc, y « 
—0.8 kpc in Fig.j^stands out, both in terms of age (5 Myr on av¬ 
erage) and mass (3 x 10® Mq). The on-going star formation in 
this region contrasts with the surrounding stellar material that has 
formed either at the tip of the bar a few ~ 30-40 Myr before, 
or even earlier elsewhere along the bar or in the spiral arms. The 
reason for this is the interaction of two gas clouds in this region. 
Fig. [TT] shows the distribution of gas, in the densest regions of the 
clouds, which exhibits tail features. These tails are a signature of 
the past and still on-going tidal interaction between the clouds, in 
a very similar fashion as interacting galaxie s exhibit long tails at 
much larger scale (e.g. IDuc & Renaudl2()I^ . 

The Eulerian nature of the code used makes it difficult to track 
the cloud progenitors back. However, by tracking back the stel¬ 
lar particles found in this structure, we are able to get insights on 
the formation of this object. Fig. [H shows the identification of the 
clouds in earlier snapshots, respectively 30 and 15 Myr before the 
instant of Fig.j^and Fig.[TT] The two progenitors visible in Fig. m 
are labelled “PI” and “P2” here. (Distinguishing which is which on 
Fig.dXjis impossible because of a too low time sampling of the sim¬ 
ulation snapshots compared to their orbital timescale around each 
other, and because their stellar contents are too spatially extended, 
such that they overlap on the late snapshot.) 

P2 is itself the result of a merger of (at least two, but possi¬ 
bly more) beads-on-a-string clouds that formed along a spiral arm. 
By visual identification with the real Galaxy, these clouds would 
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Figure 12. Positions of the cloud progenitors of the cloud-cloud collision, 30 Myr (left) and 15 Myr (right) before the instant showed in Fig. \m The “PI” 
object is formed along the edge of the bar, while the other one (“P2”) originates from a beads-on-a-string structure in one spiral arm. The two meet at the 
extremity of the bar. The progenitor “P2” is itself the result of the collision of smaller clouds that occurs between the two moments pictured here. The Scutum- 
Crux and Norma spiral arms are visually identified with those of the real Galaxy, only for the sake of clarity. The label “bar edge” points to elongated gas 
structures associated to the bar, i.e. not spirals. 


be found in the Norma arm, as labelled on Fig.[T2 Approximately 
when the P2 structure crosses the har major axis, at its extremity, or¬ 
bital crowding makes its two components interact and merge. A few 
percent of the stellar content of the final object (Fig. |9j is formed 
during this event, but a significant fraction of the stellar produc¬ 
tion gets unbound and gets mixed with field stars. This is due to 
the time-varying potential, tides, asymmetric drift to a lower ex¬ 
tent, and the internal evolution of the stellar association (although 
the collisionless aspect of the code and the softening of the gravita¬ 
tional potential at the resolution of the grid introduce biases in the 
binding of the stellar structures). 

A few Myr later, the resulting P2 is found on an orbit very 
close to that of PI, another cloud that formed along the edge of 
the bar. The two progenitors are then accelerated along the leading 
side of the bar (bottom on the right-panel of Fig.ll2t with a relative 
velocity of 85 km s“^. While leaving the tip of the har, the clouds 
gravitationally capture each other and interact for the first time (t — 
795 Myr), 5 Myr before the instant of Fig. M which corresponds 
to the mean age of the stars in the resulting cluster (Fig. H. The 
combined mass of the two clouds reaches 6 x 10® Mq, i.e. the 
regime of molecular complexes. The relative velocity of the two 
proge nitors is higher th an that observed in the Milky Way spiral 
arms dFukui et al.ll2014l) . highlighting again the peculiar dynamics 
of the tips of the bar. 

Fig. [T^ shows that the interaction of the two progenitors in¬ 
duces a rapid rotation of the merger (anti-clockwise on Fig. [ID in 
their orbital plane, i.e. the plane of the galactic disc. The innermost 
~ 3 pc of the merger host a complex velocity field resulting from 
the overlap of the internal fields of the two dense structures visi¬ 
ble in Fig. dD as well as that of the tidal bridge connecting them. 
At larger distances, we find tangential velocities of the order of 
30 km (and locally up to 100 km s“^), followed by a Keple- 
rian decrease up to the point (~ 100 pc) where kpc-scale dynam¬ 
ics takes over in setting the velocity field. Conversely, a “control” 



distance to cloud centre [pc] 


Figure 13. Tangential velocity of the gas surrounding the cloud-cloud 
merger shown in Fig. [n compared to that of a “control” cloud which has 
not experienced any recent interaction. Distances and velocities are com¬ 
puted in the reference frame of the cloud centre of mass. The merger ex¬ 
hibits a rapid rotation (with a complex behaviour in the central ~ 3 pc), 
induced by the cloud-cloud interaction, while the rotation of the control 
cloud is several times slower. 

cloud of comparable central density, in the same environment, and 
which has not experienced any recent interaction, does not exhibit 
such rapid rotation. Instead, motions around the control cloud are 
mainly radial (~ 20 km s“^), reflecting the on-going accretion of 
surrounding gas. 

During the cloud-cloud collision event, the increase in density 
and pressure triggers the bulk of star formation over ~ 20 pc, while 
the densest parts of the clouds continue spiralling around each 
other (as visible on Fig. ED- This agrees with the observational and 
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Figure 14. Schmidt-Kennicutt diagram of all the clouds in the simulation, 
highlighting the progenitors and their merger of the cloud-cloud collision 
discussed in the text. Surface densities are computed using the dense phase 
only (> 2000 cm“®) at the scale of 50 pc. The dashes lines indicate 
power-law fits to all the clouds, and the solid lines show the average val¬ 
ues. 


theoretical works which showed that such events could trigger an 
episode of star formation (among others, see Loren|l976yTmj200d; 
Ipurukawa et al.l 1 20091 : iTasker & Tarj l2009l : Inoue & Fukuill201^ . 
On top of increasing the star formation rate, collisions are also sus¬ 
pected to g ather the physical conditions for the formation of mas- 
sive stars dAnathpindikal l20ld : IXakahira et al.l |2014| : iMotte et"^ 
l2014h . Although this aspect cannot be probed in our simulation, 
it suggests that the kpc-scale dynamics of the bar, through cloud- 
cloud collisions, could lead to the formation of massive stars in 
excess with respect to those formed in other, more quiescent envi¬ 
ronments. This excess could either reflect a size-of-sample effect in 
which the production of more stars overall is likely to sample out 
further into the high-mass tail of the initial mass function (IMF), 
or it could reflect a real change in the IMF toward a statistically 
significant increase in the ratio of high to low mass stars. 

Fig.m shows a Schmidt-Kennicutt diagram of the pre- and 
post-collision stages, comparing the progenitors and the merger to 
all the clouds in the simulation. By applying a density selectiorQ 
above 2000 cm“®, we find that, before the collision, PI is among 
the most massive clouds (2 x 10® Mq), while P2 is closer to aver¬ 
age (4 X 10® M 0 ). Their Mach numbers are comparable with that 
of the bulk of the clouds, i.e. in the transonic regime (~ 0.8 at the 
scale of 24 pc). After the P1-P2 collision, most of the clouds in the 
simulation have fragmented further and formed their star clusters. 
This translates into a mild shift of the Schmidt-Kennicutt relation 


^ Note that we use this approach to minimize the numerical noise and the 
effect of star formation stochasticity, which become significant at the scale 
of individual clouds in our simulation. By using a density threshold and thus 
focussing on the densest regions of molecular clouds (~ 5-10 pc), our 
method differs from observational methods where a beam of fixed size is 
used, which explains slight differences in the normalisa t ion of the relations 
shown in Fig. 1141 and the literature. See iKraliic et alj j2014h for a more 
observation-oriented approach. 


toward higher surface densities of gas and SFR. The remnant of the 
P1-P2 collision weights 5 x 10® Mq (in dense gas), and yields a 
higher Mach number than its progenitors (« 1.3). This translates in 
the Schmidt-Kennicutt diagram as the merger hosting a more effi¬ 
cient star formation (i.e. a shorter gas depletion t ime) t han average. 
This confirms the predictions of [Renaud et alj ( 1201 2h that super¬ 
sonic turbulence increases the efficiency of the conversion of gas 
into stars or in other w ords, decreases the gas depletion time (see 
also IKraliic et alj[2014h . We note that such evolution is compara¬ 
ble to that of galaxies jumping from the regime of is olated discs to 
that o f starbur st during a galaxy -galaxy interaction teenaud et al.l 
|2014 see also lOaddi et al. l201oiff We also note that the morphol¬ 
ogy of collapsing beads-on-a-string seen in the spiral arms in our 
simulation suggests that most star formation over the galactic disc 
is spontaneous and not triggered by a collision between two large 
clouds. 

In term of complex formation history, extreme mass, and be¬ 
cause of its location close to the tip of the bar, our object re¬ 
sembles the observed complex W43 which is thou ght to origi¬ 
nate fr om multiple cloud accretion events. In particular. lMotte et al.l 
(I 2 OI 4 I) proposed that W43 could have been formed through the ac- 
cretion of clouds for m the nearby Scutum-Centaurus spiral arm. 
ICarlhoff et al.l (1201 3h noted a strong rotation of the main cloud of 
the W43 complex, together with elongated, arm-like, structures. 
These are further hints of a recent collision, in line with the forma¬ 
tion scenario proposed from our simulation. We note that the orbit 
of our W43-analogue maintains it along the edge of the bar and at 
the tips, so that it does not experience the st rong shear seen inside 
the bar. This confirms the suggestion from ICarlhoff et akl ( l2013l) 
that W43 is stable against shear, and could survive several 10 Myr 
before being destroyed by stellar feedback or tidal effects. 

Despite a few difference|3> the formation of our object 
yields remarkable similarities to that of W43. In that sense, with¬ 
out providing a perfect match, our simulation would picture a 
W43-analogue, confirming the possibility of the scenario invok¬ 
ing accretion and collision of (at least) two clouds. In our sim¬ 
ulation, some events taking place before the main collision it¬ 
self have led to the formation and ejection of several loosely 
bound or unbound stellar associations over the last 30 Myr. 
A fraction of this stellar population might be identified as the 
red super giant clusters (RSGC) detected in the vicinity of W43 
dNegueruela et alJl2010L 2012l:lGonzalez-Fernandez & Negueruelal 


I 2 OI 2 I : iGonzalez-Femandez^alj 20121) . This would then confirm 
the hypothesis of Motte et alj ( 12014 ) that these observed objects 
(« 20 Myr-old) originate from previous star forming events in 
W43 and have decoupled from the molecular complex at an early 
stage. A more detailed analysis of this point would require a colli- 
sional treatment of gravitation, which is not included in our simu¬ 
lation. 

Although cloud-cloud interaction are frequent at the tip of the 
bar due to orbital crowding, the formation of an object as mas¬ 
sive as W43 remains a rare event. In particular, we have not found 
a comparable structure at the other extremity of the bar despite a 


® Pushing the analogy further would require to check that the turbulence of 
the merger is dominated by its compressive mode. However, the resolution 
of our simulation forbids such analysis at the scale of molecular clouds. 

^ In particular, the observed object is younger (1 Myr) than the simulated 
one, and found closer to the extremity of the bar. 
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Figure 15. Map of the decoupling force between the stars and the gas 
clouds, overlaid on the gas density map (in greyscale) at f = 800 Myr. 
Red areas indicate that stars are efficiently decoupled from the cloud in 
which they formed. (See text and Appendix IaI) 


strong symmetry of the large scale galactic features (bar and spiral 
arms). The same situation occurs in observations, where no W43- 
equivalent has been detected at the remote extremity of the bar. 


3.5 Cloud-star decoupling 

3.5.1 Process 

iBash et alJlll977h an dH ashi 1 I 1979 I) found that stars remain associ¬ 
ated to the cloud in which they formed for 30-40 Myr. After this, 
they decouple from the gas structures, either because the cloud gets 
destroyed by feedback, or because large-scale dynamics makes the 
orbits of the stellar and gas component diverge. Such effect is in¬ 
deed present in our simula tion, as an asymme tric drift of young 
stars in the spiral arms (see [Renaud et alj|2013l . their Section 4.5). 
In the vicinity of the bar, i.e. at smaller galactic radii, the time-scale 
of the asymmetric drift is shorter, and other dynamical aspects like 
shear and tides must also be taken into account in the decoupling 
process. 

To understand to what extent young stars can be used to probe 
the physical conditions of their formation sites, we measure the 
intensity of a “decoupling force” that we define as the velocity of a 
young star (< 10 Myr) with respect to its closest cloud, divided by 
the age of the star and multiplied by its mass (see Appendix IaI for 
details). Fig.llSIshows that the decoupling force is the strongest at 
the tips of the bar and, to a lower extent, stronger along the edge of 
the bar than in spiral arms, as expected. 

On top of the usual decoupling effect due to the dissipative 
nature of the gas, kpc-scale dynamics accelerates the drift of stars 
away from their birth places. Along the edge of the bar (parallel 
to the major axis), the small galactocentric distance implies an en¬ 
hanced effect of tides. However, the higher density of (massive) 
clouds and star clusters due to orbital crowding at the extremities 
of the bar also induces a significant cloud-cloud tidal effect (with 
signatures as those shown in Fig. [TT] in the most extreme cases). 



Figure 16. Gas density at the position of the stars forming at f = 790 Myr 
(filled circles), and stars forming at i = 770 Myr (empty circles), taken at 
the time of their formation on the horizontal axis, and at t = 800 Myr on 
the vertical axis. Dotted lines indicate equality of the two densities. Only a 
fraction of stars is shown, for the sake of clarity. The top-left panel shows 
stars independently of their location in the galaxy, while the other three 
panels are subsets in specific formation regions (identified visually). 

The tidal forc^ is locally up to six times stronger in between the 
nearby clouds at the tips of the bar than along the bar edge. 

Shear in the dense regions is of comparable amplitude at the 
extremities and along the edge of the bar, but varies significantly 
when considering less dense media (< 500 cm~®), as noted by 
lEmsellem et al.l fioisl) . 

Finally, the combination of locally strong tides, shear and the 
dissipative nature of the gas enhances significantly the decoupling 
of stars from their cloud progenitors at the tips of the bar, with 
respect to any other region in the galaxy. Young stars are thus likely 
to be observed far from their gas nursery once they have passed by 
the extremities of the bar. This aspect makes the interpretation of 
observational data in this environment even more challenging than 
previously thought. 

3.5.2 Results 

In Fig. [T^ we compare the gas densities at the positions of stars 
(hereafter the ambient density) at the time of their formation with 
that 10 and 30 Myr later. In practice, we select stars younger than 
1 Myr att = 770 Myr and identify them at f = 800 Myr. They 
constitute our 30 Myr-old sample. The ambient density is that of 
the gas cell in which the star particle is found. We repeat the exer¬ 
cise by selecting young stars att = 790 Myr and identifying them 
at f = 800 Myr to get our 10 Myr-old sample. 

The evolution of the galaxy between t = 770 Myr and t = 

® computed as first order finite differences of the gravitational force at the 
scale of 24 pc. 
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Figure 17. Cumulative fraction of stars as a function of the ratio of ambi¬ 
ent gas densities between the time of star formation and 10 Myr (top) or 
30 Myr (bottom). In all cases, the majority of the stai's experience a de¬ 
crease of the ambient gas density, and this already after only 10 Myr. 


790 Myr and the on-going collapse of clouds already noticed in 
Fig.ll4ltranslate here into a rise of the average ambient density (i.e. 
an offset of the filled symbols to the right on Fig. [Ill • Fig.llVIshows 
the cumulative fraction of the stars as a function of the ambient 
density ratio (density at formation divided by density 10 or 30 Myr 
later). 

The gas density at the time of formation is fairly independent 
of the environment (spiral vs. bar). For about one third of the stars, 
the ambient density increases with time: either the star remains in 
the central region of its cloud which keeps collapsing, or the star or¬ 
bit makes it fly through the central region. For the other two thirds, 
the ambient density drops as the stars age, as a result of feedback 
dissolving the cloud, and/or of the cloud-star decoupling. As dis¬ 
cussed above, the decoupling effect is environment-dependent, and 
so is the ambient density drop. 

The picture is made more complicated by the environmental 
dependence of the density around the cloud itself. Since the scale of 
molecular clouds is comparable to the tangential extent of a gaseous 
spir al arm, once a star l eaves a cloud there, it also leaves the arm 
(see lRenaud et al.ll2013L their Fig. 14). Therefore, its ambient den¬ 
sity drops from that of the cloud centre (~ 10“*“® cm“® in our sim¬ 
ulation) to that of the inter-arm region (~ 1-10 cm“®) within a 
few 10 Myr. At the tip of the bar however, the inter-cloud medium 
remains relatively dense (-^ 10^ cm~®, recall Fig.O. Therefore, 
when a star leaves its cloud there, it experiences a less abrupt 
change of its ambient density. As a result, the stronger decoupling 
effect in the bar than in spirals (Fig. O is partially balanced by a 
less abrupt variation of the ambient density. 

Stars formed at the edge of the bar experience an intermediate 
effect between those in the bar extremities and in the spirals. The 
rapid transit of stars along the bar edges due to the fast orbital mo¬ 
tion on a: 1 -orbits also implies that they will experience the strong 
decoupling effect found at the tips within the first 20 Myr of 
their life (recall Section|T2}. 

In conclusion: compared to stars in spiral arms, a larger fac¬ 
tion of the stars formed at the tip of the bar evolve into a low density 


medium (Fig. [nil. Despite a denser inter-cloud medium, the ambi¬ 
ent density after a given period of evolution (10 or 30 Myr) is on 
average lower in the bar than in spiral (Fig. EDi. 


3.5.3 Implications 


On top of the mere divergence in spatial distribution of the young 
stars relative to their birth places, the decoupling plays a potentially 
important role on the effects of feedback. In all regions, the decou¬ 
pling effect is already measurable in the first 10 Myr of the star life, 
i.e. even before the supernova (SN) blasts (set at 10 Myr after star 
formation in our simulation). This implies that the feedback energy 
is not injected in the densest regions of clouds, but rather at their 
periphery or even in the inter-cloud medium. SN bubble s can then 
expan d more easily in a relatively low density medium jChevaliej 
Gia), and potentially drive galactic fountains or outflows. At first 
order, the net feedback effect on the ISM thus depends mainly on 
two factors: (i) how fast a star decouples from its cloud (Fig. [Bl, 
and (ii) how dense the ISM surrounding the cloud is (Fig. Ei and 

EU. 


This decoupling aspect is not retrieved in most of the other 
simulations. On the one hand, simulations of the ISM or of isolated 
clouds miss the large-scale dynamics of the galaxy and thus the 
very drivers of the decoupling (asymmetric drift, shear, tides). On 
the other hand, the decoupling scales (-^ 10 pc) are not resolved 
by the simulations of large cosmological volumes. The injection of 
feedback is thus inaccurate in both cases. 

The decoupling effect in bars that we noted could play a role 
in removing low angular momentum gas via outflows ( or fountains) 
at ear ly stages of galaxy formation, as proposed by ISrook et'5] 
( l2012h . The low angular momentum gas, sitting at small galactic 
radii, ultimately gives rise to overly massive bulges and peaked ro¬ 
tation curves, in contrast to observed galaxies. This so-called “an¬ 
gular momentum problem” has plagued galaxy formation simula- 


a gued g ali 

tions for decades (e.g . |Abadj_eLal][20()3|; |Nayarro_^White 
Govemato et^ 2007: l^annapieco et al.ll2009l: Ubler et al 


1994 


2014 


Agertz & Kray tsovll 20151) . Strong stellar feedback provides a solu 

tion, although t his may have a det rimental effect on galactic thin 
disk formation dRoskar et al.ll2014h . By allowing for SN feedback 
to more efficiently blow out pockets of low density gas in specific 
galactic regions, cloud-star decoupling in bars could be an impor¬ 
tant ingredient in alleviating the problem. 

Reality is much more complex than the picture we draw up 
here. The ISM is clumpy and non-uniform, making feedback bub¬ 
bles asymmetric. Furthermore, clustered star formation implies that 
multiple, asynchronous, blast waves would interact in a given star 
forming volume. On top of this, colliding feedback shells and 
interaction with nearby clouds are likely to trigger star forma¬ 
tion and thus subsequent injection of feedback in an already al¬ 
tered medium. A full understanding of the physics of feedback 
in galactic (or even cosmological) context requires more simula- 
tions like our s , whic h is still very costly. The zoom-in method by 
iBonnell e t al.l (l20 13h is likel y to provide insights in this topic (see 
also lRev-Raposo et al.ll201^ . 
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4 DISCUSSION AND CONCLUSIONS 

We use a hydrodynamical simulation of a Milky Way-like galaxy to 
probe the role of the har environment in triggering and regulating 
the formation of clouds and stars. Our approach consists in con¬ 
necting kpc-scale processes with the formation and early evolution 
of parsec-scale structures. Our main conclusions are: 

• The leading edges of the har yield converging flows and large- 
scale shocks that favour the gath ering of dense gas in kpc-long fil¬ 
aments, confirming the results of lAthanassou 3 (Il992h and others. 
Furthermore, these regions host supersonic turbulence which then 
helps the fragmentation of these filaments into clouds. 

• The tangential velocity gradient at the edge of the bar is 
much weaker than in innermost regions, while the radial compo¬ 
nent yields shocks. Such differences explain that gas clumps are 
likely to survive at the edge but not inside the bar. 

• The orbital pattern in the har creates a fast circulation of matter 
along the har. While going from one tip to the other takes only 
20 Myr, the gas slows down at both tips and accumulates there. 

• Because of the high velocities along the bar, the objects older 
than 10 Myr have formed in a different environment and under 
different physical conditions than the ones they are when observed. 

• Orbital crowding at the tip of the bar leads to cloud-cloud col¬ 
lisions that can form massive molecular complexes like W43 and 
massive stellar associations (~ 10® Mq). Such objects remain at 
the tip of the bar for a rather short period (~ 10 Myr or less) after 
their formation. 

• A cloud-cloud collision increases the rotation of the merger, 
the turbulent Mach number inside the cloud and the efficiency of 
star formation, in a comparable way as galaxy mergers. 

• Young stars decouple from their gas cloud more efficiently at 
the tips of the bar than along the edge or in spiral arms. Such aspect 
can affect the net effect of stellar feedback, possibly even at galactic 
scale. 

All together, these points illustrate the paramount effect of the 
fast dynamical evolution connected to the bar: stars even only a 
few 10 Myr old can be found in a very different environment than 
the one they formed in. Tips of the bar play an important part in 
this picture. Either they host clouds long enough to witness their 
collapse and star formation, or they gather physical conditions al¬ 
lowing star formation events (like cloud-cloud collisions). In the 
latter case, the actual event can occur a few Myr later, while the 
clouds already exited this particular region. 

The pro cesses we highligh t are likely to be altered in weaker 
bars (but see Izhou et al.l 20151 and references therein) or by the 

S iresence of multiple b ars, especially if they rotate at different speed 
Wozniak et al.l[l99^ . e.g. by changing the large scale velocity 
field. Active star forming regions, comparable to those we anal¬ 
ysed, has been obse rved at the extremities of bars of s everal galax¬ 
ies, likeNGC 1097 ( fOndrechen & v an der Hul sd 19831) . NGC 1433, 
NG C 1512 and NGC 53 83 jAthanassoulJl99jk 

iKraliic et al.l (1201 2 h showed that, although bars are short-lived 
at high redshift {z > 1), those formed at a « 0.8-1 survive for 
long periods. They evolve through the accretion of intergalactic 
gas that they re-distribute with in the galactic disc via gravitational 
torques jBoumaud et al.|[2005tl . Acc ording to the classifi cation of 
bars in isolated galaxies proposed bv IVerlev et alj l l2007h . the case 
modelled here is representative of t he most common c lass (their 
group E). The evolution scenario in IVerlev et alj l l2007h suggests 


such bar would weaken and even be destroyed by gas infall. There¬ 
fore, simulations at parsec resolution and in cosmological context 
are needed to establish the frequency and the importance of the as¬ 
pects we presented here in the special case of our Milky Way-like 
simulation. 
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APPENDIX A: COMPUTATION OE THE DECOUPLING 
EORCE 

To measure the effect of the decoupling of the stars from the gas 
cloud in which they form, we first identify clouds as density peaks 
in one snapshot of the simulation. We select the grid cells in which 
the gas is denser than the star formation threshold (2000 cm“®). 
All contiguous selected cells constitute one cloud. Then, we select 
stars between 2 and 10 Myr old, and identify their closest cloud 
as their formation site. Choosing older stars could lead to errors in 
the identification of the formation site. We compute the velocity of 
the star with respect to the cloud and its radial and tangential com¬ 
ponent. We reject the stars with a negative radial velocity, as they 
are moving toward the cloud and do not experience a decoupling 
effect. The total velocity is then converted into a force, by multi¬ 
plying it by the mass of the star, and dividing it by the age of the 
star. Einally, we compute the average force for each cloud over all 
its stars. 

To obtain the map of Pig. [B] we convolve the scattered dis¬ 
tribution of the force for each cloud by a Gaussian kernel. We then 
multiply this map by a smoothed version of the gas density map for 
graphical purposes only. 
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